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We explore leptogenesis induced by the propagation of neutrinos in gravitational backgrounds that 
may occur in string theory. The first background is due to linear dilatons and the associated Kalb- 
Ramond field (axion) in four non-compact space-time dimensions of the string, and can be described 
within the framework of local effective lagrangians. The axion is linear in the time coordinate of 
the Einstein frame and gives rise to a constant torsion which couples to the fermion spin through a 
gravitational covariant derivative. This leads to different energy-momentum dispersion relations for 
fermions and anti-fermions and hence leptogenesis. The next two backgrounds go beyond the local 
effective lagrangian framework. One is a stochastic (Lorentz Violating) Finsler metric which again 
leads to different dispersion relations between fermions and antifermions. The third background 
of primary interest is the one due to populations of stochastically fluctuating point-like space-time 
defects that can be encountered in string/brane theory (DO branes). Only neutral matter interacts 
non-trivially with these defects, as a consequence of charge conservation. Hence, such a background 
singles out neutrinos among the Standard Model excitations as the ones interacting predominantly 
with the defects. The back-reaction of the defects on the space-time due to their interaction with 
neutral matter results in stochastic Finsler- like metrics (similar to our second background). On 
average, the stochastic fluctuations of the DO brane defects preserve Lorentz symmetry, but their 
variance is non-zero. Interestingly, the particle/antiparticle asymmetry comes out naturally to 
favour matter over antimatter in this third background, once the effects of the kinematics of the 
scattering of the D branes with matter is incorporated. 
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I. INTRODUCTION 

One of the most important issues of fundamental physics, relates to an understanding of the magnitude of the 
observed baryon asymmetry rig — n-g (where B denotes baryon, B denotes antibaryon, ng is the number density 
of baryons and n-g the number density of antibaryons in the universe). The numbers of protons and neutrons 
far exceeds the number of antiprotons and antineutrons.The universe is overwhelmingly made up of matter rather 
than anti-matter. According to the standard Big Bang theory, matter and antimatter have been created in equal 
amounts in the early universe. However, the observed charge-parity (CP) violation in particle physics prompted 
A. Sakharov [2] to conjecture that for Baryon Asymmetry in the universe (BAU) we need: 

• Baryon number violation to allow for states with AB ^ starting from states with AB = where AB is the 
change in baryon number. 

• If C or CP conjugate processes to a scattering process were allowed with the same amplitude then baryon 
asymmetry would disappear. Hence C and CP need to be broken. 

• Chemical equilibrium does not permit asymmetries. Hence Sakharov required that chemical equilibrium does 
not hold during an epoch in the early universe. 

Hence non-equilibrium physics in the early universe together with baryon number (B), charge (C) and charge-parity 
(CP) violating interactions/decays of anti-particles in the early universe, may result in the observed BAU. In fact 
there are two types of non-equilibrium processes in the early universe that can produce this asymmetry: the first type 
concerns processes generating asymmetries between leptons and antileptons (leptogenesis) , while the second produces 
asymmetries between baryons and antibaryons (baryogenesis) . The near complete observed asymmetry today, is 
estimated in the Big-Bang theory |3j to imply: 

An(T ~ 1 GeV) = nB ~ n B „ Z n s = (8 . 4 _ 8 9) x 10 -u (1) 

at the early stages of the expansion, e.g. for times t < 10~ 6 s and temperatures T > 1 GeV. In the above formula Ub 
(n-g) denotes the (anti) baryon density in the universe, and s is the entropy density. Unfortunately, the observed CP 
violation within the Standard Model (SM) of particle physics (found to be of order e = O(10 -3 ) in the neutral Kaon 
experiments [l|) induces an asymmetry much less than that in (JXJ) [4|. There are several ideas that go beyond the SM 
(e.g. grand unified theories, supersymmetry, extra dimensional models etc.) which involve the decays of right handed 
sterile neutrinos. For relevant important works on this see [BI-ITTj . These ideas lead to extra sources for CP violation 
that could generate the observed BAU. Fine tuning and somewhat ad hoc assumptions are involved in such scenarios 
and the quest for an understanding of the observed BAU still needs further investigation. 

The requirement of non-equilibrium is on less firm ground than the other two requirements of Sakharov, e.g. if 
the non-equilibrium epoch occurred prior to inflation then its effects would be hugely diluted by inflation. A basic 
assumption in the scenario of Sakharov is that CPT symmetry [l2T ] (where T denotes time reversal operation) holds in 
the very early universe which leads to the production of matter and antimatter in equal amounts. Such CPT invariance 
is a cornerstone of all known local effective relativistic field theories without gravity, which current particle-physics 
phenomenology is based upon. It should be noted that the necessity of non-equilibrium processes in CPT invariant 
theories can be dropped if the requirement of CPT is relaxed. This violation of CPT (denoted by CPTV) is the result 
of a breakdown of Lorentz symmetry (which might happen at ultrahigh energies [13]). For many models with CPTV, 
in the time line of the expanding universe, CPTV generates first lepton asymmetries (leptogenesis); subsequently 
through sphaleron processes or Baryon-Lepton (B-L) number conserving processes in grand unified theories, the 
lepton asymmetry can be communicated to the baryon sector to produce the observed BAU. 

Thus, CPTV in the early universe may also obviate the need for including extra sources of CP violation, such as 
sterile neutrinos and/or supersymmetry, in order to obtain the observed BAU. In this article we will consider CPTV 
leptogenesis from a non-Riemannian point of view, inspired by a stringy model of gravitational defects and backgrounds 
interacting with neutral fermions. The explicit non-Riemannian structure that we investigate is Finsler geometry 
[Tij where momentum as well as position are explicitly involved in its metric. The defects that we will consider 
are point-like solitonic structures in some string theories; they are known as DO-branes (or D-particles) [15]. Our 
generic model involves effectively three-space-dimensional brane uuniverseni verses, obtained from compactification of 
higher-dimensional branes, which are embedded in a bulk space punctured by D-particles (see e.g. [16]). One of 
these brane universes constitutes our observable world, which moves in the bulk. As a consequence of this motion, 
D-particles cross the brane. Open strings for electrically neutral particles on the brane can attach an end to the 
D-particle; subsequently this string can detach (l7l. Il8j. Scattering off a population of D-particles (D-foam) affects 
the kinematics of the stringy matter and leaves an imprint on the background geometry [19j on the brane world. This 
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geometry is similar to Finsler metrics but with stochastic parameters [20-22]. We first investigate the consequences 
for gravitational leptogenesis of this metric structure in the general setting of Finsler geometry without reference 
explicitly to D-foam . For the underlying stringy model however, since it is microscopic, we can consider in addition 
the kinematics of D particle scattering. This kinematical aspect leads naturally to the asymmetry between the particle 
and anti-particle abundances having the right sign. The kinematical argument involves recoil kinetic energy does not 
fit into an effective local field theory approach and represents a new approach that is relevant to leptogenesis. 

The structure of the article will be as follows: in the next section [II] we shall briefly review some relevant existing 
models for fermionic asymmetry, which entail CPTV in the early universe as alternatives to scenarios involving sterile 
neutrinos and/or supersymmetry. The discussion will also help to differentiate between our viewpoint and those 
prevailing in the literature. The models for leptogenesis that we will discuss are of gravitational type, namely CPTV- 
induced differences in the dispersion relations between particles and antiparticles propagating in these backgrounds. In 
section Hill we propose a new model for gravitational leptogenesis which follows broadly an earlier framework [23ll24j. 
but differs crucially in that the full gravitational multiplet [25| that arises in string theory is used. The leptogenesis 
in this model is due to CPTV dispersion relations between fermions/antifermions, induced by the (constant) torsion 
associated with the antisymmetric Kalb-Ramond tensor field in the gravitational multiplet. The torsion is space-time 
independent for a particular solution of the conformal invariance equations of the string associated with a linear 
dilaton. In the proposals discussed in section [TT] the space-time independence of the torsion was not guaranteed. In 
section ITVl we consider another scenario for gravitational leptogenesis that involves a non-Riemannian Finsler metric, 
with stochastically fluctuating parameters, a variant of which appears in our string/brane model in section [Vj In 
section [V] we present our model for the stringy universe and the induced CPTV and leptogenesis/baryogenesis. The 
basic properties of the model are reviewed briefly in subsection IV Al where our observable universe is represented as 
a brane world, punctured by fluctuating "point-like" brane Defects (D-particles or DO-branes); the back- reaction of 
D-particles, during their topologically non-trivial scattering with neutrinos, results in local metric distortions. These 
distortions depend, however, on both the coordinate and momentum variables of the particles, thus inducing a kind 
of Finsler geometry. In subsection IV Bl we discuss the kinematics of scattering of the D-particle (DO— brane) and 
stringy matter in addition to an induced Finsler-like metric. This together with the stochasticity of the recoil process 
leads to CPTV and matter dominance over antimatter naturally, without the need for an adjustment of the sign 
of the lepton/antilepton asymmetry. This is an interesting feature of our model which differentiates it from earlier 
proposals on gravitational leptogenesis/baryogenesis (23l. [24. [2^ - [32l | . where the sign of the asymmetry is implicitly 
chosen. In common with earlier discussions of CPTV in baryogenesis [33], this class of models involves violation of 
Lorentz symmetry, but only because of a non-zero variance of the stochastic parameter. Estimates on the range 
of parameters of the model, for which the induced asymmetry can lead to the observed Baryon Asymmetry in the 
universe, are given. We conclude that D-particles of Planck-size mass are probably needed for this purpose, if standard 
cosmological properties of neutrinos are assumed. Conclusions and outlook are finally given in section IVT1 



II. MODELS OF CPT VIOLATION (CPTV) 

In this section we shall review some existing models of CPTV induced asymmetry between matter and antimatter 
in the early universe, which can be contrasted with our approach in this article. We shall be brief in our exposition, 
referring the interested reader to the relevant literature for more details. 



A. CPTV Models with Particle-Antiparticle Mass Difference 

The simplest possibility [34] for inducing CPTV in the early universe is through particle-antiparticle mass differences 
to to. These would affect the particle phase-space distribution function /(i?,/i) 

/( J B, M ) = [exp( J B-M)/T)±l]- 1 , E 2 = f + to 2 (2) 

and antiparticle phase-space distribution function 

i ^2 -e , -2 



f(E,fi) = [exp(E-fi)/T)±l}- 1 , E J =0 



m 



with p being the 3— momentum. (Our convention will be that an overline over a quantity will refer to an antiparticle, + 
will denote a fermionic (anti-)particle and — will denote a bosonic (anti-)particle.) Mass differences between particles 
and antiparticles, to — m ^ 0, generate a matter-antimatter asymmetry in the relevant densities 

n-n = g d .o.f. J [f(E, ft) - f(E, p)] , (3) 
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where gd.o.f. denotes the number of degrees of freedom of the particle species under study. In the case of spontaneous 
Lorentz violation [33j there is a vector field with a non-zero time-like expectation value which couples to a global 
current J M such as baryon number through an interaction lagrangian density 

C = XA^. (4) 

This leads to m ^ fh and /i 7^ p,. Alternatively, following [34j we can make the assumption that the dominant 
contributions to baryon asymmetry come from quark-antiquark mass differences, and that their masses "run" with 
the temperature i.e. m ~ gT (with g the QCD coupling constant). One can provide estimates for the induced baryon 
asymmetry on noting that the maximum quark-antiquark mass difference is bounded by the current experimental 
bound on the proton-antiproton mass difference, 5m p (= \m p — m p |), known to be less than 2 • 10 -9 GeV. Taking 
rt 7 ~ 0.24 T 3 (the photon equilibrium density at temperature T) we have |34| : 



Ar = — = 8.4 x 10 — , Sm q = \m q - m q \ . (5) 

is too small compared to the observed one. To reproduce the observed (3t=o ~ 6 • 10~ 10 one would need 8m q (T = 
100 GeV) ~ 10~ 5 — 10~ 6 GeV > Sm p , which is somewhat unnatural. 

However, active (light) neutrino-antineutrino mass differences alone may reproduce BAU; some phenomenological 
models in this direction have been discussed in (35j . considering, for instance, particle-antiparticle mass differences 
for active neutrinos compatible with current oscillation data. This leads to the result 

n B =n v -n-~ — - — (6) 

yielding ns/s ~ ~ 10 -11 at T ~ 100 GeV, in agreement with the observed BAU. (Here s, n v , and/i^ are the 
entropy density, neutrino density and chemical potential respectively.) 



B. CPTV Decoherence Models 

But particle-antiparticle mass differences may not be the only way by which CPT is violated. As discussed in [36l . [37| . 
quantum gravity fluctuations in the structure of space-time, may be strong in the early universe; the fluctuations may 
act as an environment inducing decoherence for the (anti-)neutrinos. However the couplings between the particles and 
the environment are different for the neutrino and antineutrino sectors. Once there is decoherence for an observer with 
an energy for a low-energy (compared to the Planck scale Mp ~ 10 19 GeV) observer, the effective CPT symmetry 
generator may be ill-defined as a quantum mechanical operator, according to a theorem by R. Wald (38| . leading to 
an intrinsic violation of CPT symmetry. This type of violation may characterise models of quantum gravity with 
stochastic space-time fluctuations due, for instance, to gravitational space-time defects, as is the case of certain 
brane models Qil [III H2 ■ I n such a case, a slight mismatch in the strength of the stochastic space-time fluctuations 
between particle and antiparticle sectors, can lead to different decoherence parameters to describe the interaction of 
the gravitational environment with matter. 

In [36, 37], simple models of Lindblad decoherence [3^|, conjectured to characterise quantum-gravity-induced CPTV 
decoherent situations 0, |4l|, have been considered for neutrinos f42| . It was assumed on phenomenological grounds, 
that non-trivial decoherence parameters were only present in the antiparticle sector: this is consistent with the lack of 
any experimental evidence to date [43j-|45| for vacuum decoherence in the particle sector. The antineutrino decoherence 
parameters (with dimension of energy) had a mixed energy dependence. A diagonal Lindblad decoherence matrix for 
three-generation neutrinos requires eight coefficients Some of the eight coefficients were assumed for simplicity in 
[H, [33| to be proportional to the antineutrino energy 

li = in-E, { = 1,2,4,5 
11 M P 

while the remaining (subdominant) ones were inversely proportional to it 

_ 10- 24 (GeV) 2 . 

7j = ^ > J = 3, 6, 7, 8 . 

The model was proposed without any microscopic justification; its choice was originally motivated by fitting the 
LSND "anomalous data" in the antineutrino sector [46] with the rest of the neutrino data, and this required T to be 
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T/Mp ~ 10 18 ,i.e. in the temperature range of electroweak symmetry breaking. One can derive [3a, [37| an active 
(light) v — V asymmetry of order 

A = ny-ny = 7i = T_ _ # ^ 
n v + n- y^W 2 " M P ' y/Am? ' 

where Am 2 denotes the (atmospheric) neutrino mass squared difference, which plays the ri^cele of a characteristic low 
mass scale in the problem. This lepton number violation is communicated to the baryon sector by means of baryon 
number (B) plus lepton number (L) conserving sphaleron processes. These processes lead to an estimate [Hj] for the 
current value of B to be 

s g* n 7 

with n 7 the photon number density, g* the effective number of degrees of freedom (at the temperature where the 
asymmetry developed, i.e. the electroweak symmetry breaking temperature in the model of [36|). g* depends on the 
matter content of the model (with a typical range g* € [10 2 — 10 3 ]). For such parameter values A ~ 10~ 6 and so 
the observed BAU may be reproduced in this case without the need for extra sources of CP violation e.g. sterile 
neutrinos. Such models, however, do not provide an underlying microscopic understanding. In particular there is 
missing an understanding of the preferential role of the neutrino compared to other particles of the Standard Model in 
the CPT violating decoherence process. Within some microscopic models of space-time foam, involvin g po pulations 
of point-like brane defects (D-particles) puncturing three(or higher)-spatial dimension brane worlds fl6l Il8l.[2l|. such 
a preferred role may be justified as we shall discuss in section [V] Moreover D-particles also imply that the framework 
of Ricmannian geometry will need to be generalised to Finsler geometries. Finsler geometries (with a stochastic 
background) will be investigated in order to evaluate the possibilities of CPTV occuring. 



C. CPTV-induced by Curvature effects in Background Geometry 

Although the role of gravity was alluded to in the last subsection, associated features of space-time were not 
discussed. In the literature the role of gravity has been explicitly considered within a local effective action framework 
which is essentially that of (fj| A coupling to scalar curvature 1Z [HI, H3, [H, through a CP violating interaction 
Lagrangian C: 

C = -^l j d 4 x^g-(d,K)J» (9) 

where M* is a cut-off in the effective field theory and J' 1 could be the current associated with B — L ( L being the 
lepton number) . There is an implicit choice of sign in front of the interaction (0) , which has been fixed so as to ensure 
matter dominance. 

It has been shown that [27| 

s M*T d 

Td being the freeze-out temperature for B — L interactions. The idea then is that this asymmetry can be converted 
to baryon number asymmetry provided the B + L electroweak sphaleron interaction has not frozen out. To leading 
order in M~ 2 we have R — 8ttG(1 — 3w)p where p is the energy density of matter and the equation of state is p = wp 
where p is pressure. For radiation it) = 1/3 and so in the radiation dominated era of the Friedmann- Robertson- Walker 
cosmology R = 0. However w is precisely 1/3 when T£ = 0. In general T£ oc fi{g)F^ u F^ v where (3(g) is the beta 
function of the running gauge coupling oin a SU(N C gauge theory with iV c colours. This allows w =^ 1/3. Further 
issues in this approach can be found in [H [13, H, [30. 

Another approach involves an axial vector current 3, H^, HH, HH instead of J p . The scenario is based on the well 
known fact that fermions in curved space-times exhibit a coupling of their spin to the curvature of the background 
space-time. The Dirac Lagrangian density of a fermion can be re-written as: 

C = V~9^ (n a d a -rn + n 5 B a ) 1> , B d = e abcd e bx (8 a e\ + e v c e» a ) , (11) 



in a standard notation, where e M a are the vielbeins, is the Christoffel connection and Latin (Greek) letters denote 
tangent space (curved space-time) indices. The space-time curvature background has, therefore, the effect of inducing 
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an "axial" background field B„ which can be non-trivial in certain anisotropic space-time geometries, such as Bianchi- 
type cosmologies [U, [2^, [3l], [32j] . For an application to particle-antiparticle asymmetry it is necessary for this axial 
field B a to be a constant in some local frame. The existence of such a frame has not been demonstrated. As before 
if it can be arranged that B a ^ for a — then for constant Bo CPT is broken: the dispersion relation of neutrinos 
in such backgrounds differs from that of antineutrinos. Explicitly we have 



(p— B) 2 +m 2 + B a , E = \J{p + B) 2 +m 2 - B Q . (12) 

The relevant neutrino asymmetry emerges on following the same steps (c/. ([3])) used when there was an explicit 
particle-antiparticle mass difference, As a consequence the following neutrino-antineutrino density difference is found 
in Bianchi II Cosmologies (2I I2H l3lH32T| : 

An v = n v -n-^g*T 3 {^j (13) 

with g* the number of degrees of freedom for the (relativistic) neutrino. An excess of particles over antiparticles is 
predicted only when Bq > 0, which had to be assumed in the analysis of 0, HH HH, [32J; we should note, however, 
that the sign of Bq and its constancy have not been justified in this phenomenological approach Q. 

At temperatures T < Td, with Td the decoupling temperature of the lepton-number violating processes, the ratio 
of the net Lepton number AL (neutrino asymmetry) to entropy density (which scales as T 3 ) remains constant, 

AL(T <T d ) = ^~^ (14) 

which, for Td ~ 10 15 GeV and Bq ~ 10 5 GeV, implies a lepton asymmetry (leptogenesis) of order AL ~ 10~ 10 , in 
agreement with observations. The latter can then be communicated to the baryon sector to produce the observed 
BAU (baryogenesis) , either by a B-L conserving symmetry in the context of Grand Unified Theories (GUT), or via 
(B + L)-conserving sphaleron processes, as in the decoherence- induced CPT Violating case of [3(1 [33], mentioned 
previously. 

In closing this section, let us recapitulate what we regard as some of the less than satisfactory issues with the 
above proposals. The simple assumption of mass differences between particles and antiparticles dominated by the 
quark/antiquark mass differences cannot reproduce the observed BAU, because the latter are naturally bound by 
the current proton/anti-proton mass difference. Moreover, phenomenological models involving neutrino/antineutrino 
mass differences, although capable of reproducing the observed BAU, nevertheless are ad hoc and lack microscopic 
justification. In addition, the models involving spin/ curvature coupling require constancy of the axial vector B 11 , 
which is assumed without proof. The same holds for the sign of the asymmetry between matter and antimatter in 
these models. The models we will propose below address these problems in a more microscopic way and offer partial 
resolution. 



III. CPTV-INDUCED IN (STRING-INSPIRED) BACKGROUND GEOMETRY WITH TORSION 

In this section we would like to discuss an explicit case where a constant B° "axial" field arises as a consequence of 
the interaction of the fermion spin with a background geometry with torsion. This is a novel result that, to the best 
of our knowledge, has not been noticed before. In the case of torsion the Christoffel symbol contains antisymmetric 
parts in its lower indices, so that r A ^ ^ \u- As a consequence, the last term of the right-hand side of the definition 
of the B d "vector" (fTTj). which would vanish in a torsion free Einstein manifold, is not zero. Such pure torsion terms 
may then contribute to CPTV dispersion relations, as we shall now demonstrate, even if the torsion-free metric does 
not. 

The case we shall examine below is inspired from string theory. It is known that the massless gravitational multiplet 
of strings contains the dilaton (spin 0, scalar), 4>, the graviton (spin 2, symmetric tensor), <? M „ = g vil and the spin 
1 antisymmetric tensor £? M „ = —B vp fields. The (Kalb-Ramond) field B appears in the string effective action only 
through its totally antisymmetric field strength, H^ up — d^B vp ^ where [. ..] denotes antisymmetrization of the 
respective indices. It is known from string amplitude calculations [47| that H^p plays the role of (Kalb-Ramond) 



The above considerations concern the dispersion relations for any fermion, not only neutrinos. However, when one considers matter 
excitations from the vacuum, as relevant for leptogenesis, we need chiral fermions to get non trivial CPTV asymmetries in populations 
of particle and antiparticles, because < ip^'y B ip >= — < Vl7 5 V , l > + < V , Jj7 5 V'.R >• 
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torsion in a generalised connection. Indeed, the four-dimensional bosonic part of the 0(a') effective action of the 
string (in the so-called Einstein frame, where the scalar curvature term is canonically normalised) reads: 

' d i x^(R(g)~2d^d^--e-^H fwp H^ + ... ) (15) 



16?r J v "V wy p 12 

where M s — \j \fa! is the string mass scale and denotes the (dimensionless) compactification volume, where the 
compact radii are expressed in units of \fa! . From this form it becomes immediately obvious that the Kalb-Ramond 
tensor H 2 terms can be assembled together with the Einstein scalar curvature term R(g) in a generalised curvature 
R(g, r) term defined with respect to a generalised Christoffel symbol (connection) T with mixed symmetry in its lower 
indices 



pA , —2$ ttX _ pA , rp\ /-i o\ 



where r A = r „„ is the torsion- free Einstein- metric connection, and T'* = —T Vjl is the torsion. It can be 
shown that this result persists in higher orders of a', after appropriate field redefinitions, which leave the scattering 
amplitudes invariant. Thus, for our purposes below, we consider it as an exact result, valid to all orders in stringy 
cr-model perturbation theory. Fermions then in such effective theories will couple to the TJ-tensor via spin connections 
with torsion, i. e. the relevant Lagrangian terms (to lowest order in a' ) will be of the form (fTT|F|. 

In ref. [50j exact solutions to the conformal invariance conditions (to all orders in a') of the low energy effective 
action of strings have been presented. In four "large" (uncompactified) dimensions of the string, the antisymmetric 
tensor field strength can be written uniquely as 

H^p = e 2 ^^ vpa d a b(x) (17) 

with eoi23 = \fg and e^ w P a = \g\~ 1 e tiupcn with g the metric determinant. The field b(x) is a "pseudoscalar " axion-like 
field. It is worthy of mentioning that both the dilation $ and axion b fields are fields that appear as Goldstone 
bosons of spontaneously broken scale symmetries of the string vacua, and as such are exactly massless classically. In 
the effective string action such fields appear only through their derivatives hence a solution that is linear in time for 
these two fields will only shift the various minima of all other fields in the effective action that couple to them by a 
space-time independent amount. 

The exact solution of [50] is precisely such that in the string frame both dilation and axion fields are linear in target 
time X°, Q(X°) ~ X°, b(X°) ~ X°. In the "physical" Einstein frame (where cosmological observations are made) the 
temporal components of the metric are normalised to goo = +1 by an appropriate change of the time coordinate. In 
this setting, the solution of [50] leads to a linearly expanding Friedmann-Robertson- Walker (FRW) metric, with scale 
factor a(t) ~ t, with t the FRW cosmic time. Moreover, the dilaton field $ behaves as — \nt + <fro, with <fio a constant, 
and the axion field b(x) is linear in time t, that is: 

b(x) = V2e" l '« y/Q 2 ^t , (18) 

where M s is the string mass scale, n is a positive integer, associated with the level of the Kac-Moody algebra of 
the underlying world-sheet conformal field theory model and Q 2 > is its central-charge deficit (supercritical string 
theory), with the central charge being given by: c = 4— 12Q 2 — +cj, where cj is the central charge associated with 
the world-sheet conformal field theory of the compact "internal" dimensions of the string model [50] . The requirement 
of cancellation of the world-sheet ghosts that appear as a result of fixing reparametrisation invariance of the world- 
sheet coordinates forces the constraint c = 26. In the presence of a non-zero Q 2 there is an additional dark energy term 
in (fl~5|) of the form J d A x^J~ ge 2 *(— AQ 2 )/a' . It should be noted that the linear axion field (fT8|) remains a non-trivial 
solution of the string-effective-action equations even in the static space-time limit with a constant dilaton field [50] . In 
such a case the space time is an Einstein universe with positive cosmological constant and constant positive curvature 
proportional to 6/(n + 2). 

In the above solutions, the covariant torsion tensor e _2 *H All/p is constant, as becomes evident from (|16[) . (j!7[) since 
the exponential dilation factors cancel out in the relevant expressions. Only the spatial components of the torsion are 
non zero in this case, 

T ijk - e ijk b = e^k^Q^e-* -± , (19) 

\ n 



2 We note that fermions coupled to Kalb-Ramond torsion tensors H^ vp have been considered in the literature |48| , |49H but from a different 
perspective than ours. 
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where the overdot denotes time t derivative. From (flU)) . ([T7]) and (jTTj) . we observe in this case that only the temporal 
component of the B d vector has non trivial, B° ^ 0, and it is the -torsion parts that contribute to it. The torsion- 
free gravitational parts on the other hand (for the FRW or flat case) yield vanishing contributions. From (fTTj) and 
(UHJ) then we obtain a constant B° of order 

B° ~ y^e"^ ^ GeV > 0. (20) 

Notice here that the sign of B° is fixed by string theory conventions. From the previous discussion, to get the 
phenomenologically acceptable leptogenesis in such a toy model it seems to require B° ~ 10 5 GeV. Since in string 
theory M s is a free parameter, restricted by phenomenological considerations to be higher than O(10 4 ) GeV, we thus see 
that we do not get a serious constraint on the Kac-Moody level from the requirement of CPTV-induced leptogenesis. 
Since the string coupling g s — e~*° < 1 is assumed weak for phenomenological purposes (g 2 /4:iT ~ 1/20), then 
we observe from (f2"U)) that it is mainly the central charge deficit Q 2 of the underlying conformal field theory that 
determines the order of lepton asymmetries in this model. 



IV. CPTV IN STOCHASTIC FINSLER GEOMETRIES 



Although all the models displaying CPTV that we have considered so far are based on local effective field theories, 
there is no compelling reason for a restriction to such a framework. In fact a microscopic model involving space-time 
defects based on string-brane theory suggests the use of a non-Riemannian metric background similar to that which 
occurs in Finsler geometry [Til l5lj.(This model will be discussed in a subsequent section.) Independently there has 
been much interest in Finsler geometry [52-6l| for characterising the Early universe [62l - [6(| and for descriptions of 
modified dispersion relations for particle probes pill. 1(331 l67l |68| . Finsler geometry has a metric which, in addition to 
space-time coordinates, depends also on "velocities". Lorentz symmetry is broken through some fixed vectors in the 
metric. We explore the consequences of making such vectors having components which are stochastic with possibly 
zero mean. This is a feature that arises in the defect model of D-foam pj| [U, 13 that we will consider in the next 
section. The defects stochastically fluctuate, due to both statistical and quantum stringy effects in large populations of 
such D-particles that can populate eras of the early universe. As we shall discuss below, the result of the interaction of 
neutrinos with these defects, leads to stochastically fluctuating Finsler-like metrics. However we wish to consider the 
consequences for CPTV and matter-antimatter asymmetry of this stochasticity in a general context. This underlying 
model provides our main motivation to study this class of space-times in this section and to contrast our findings on 
the induced CPTV for such cases with the corresponding ones for the D-foam model. For clarity, we commence our 
discussion with a brief reminder of the definition and properties of Finsler geometries, and then we proceed to discuss 
CPTV issues in a particular, but representative, Finsler-like geometry, which however is stochastically fluctuating. 

A Finsler geometry on a manifold M is defined in terms of a Finsler norm F (x, y), a real function of two arguments 
x and y, where x € M and y € T X M (the tangent space at x) [3, [EH- The norm F (x, y) satisfies: 

. F(x,y)^0ify^0 

• F (x, Ay) = |A| F (x, y)for Ael. 

The Finsler geometry is defined in terms of a metric g^ (x, y) which is given in terms of the Finsler norm 

, , ld 2 F 2 (x,y) 

The inverse of g^ (x, y) is represented by g^ v (x, y). It is necessary to define a similar structure in phase space (i.e. 
for the co-tangent space). The dual (y) to y M is defined by [5l| 

w m (v) = 9n» (x,y)y u - (22) 

Furthermore this relation can be inverted so that corresponding to u> there is a dual y (cj) . The Finsler norm G {x, w)on 
cotangent space is defined by 

G(x,Lj) = F(x,y(uj)). (23) 
The Finsler metric (x,u>) in phase space is defined analogously to the usual Finsler metric in (f2"Tj) 

l d 2 G 2 (x,u) 

h flOto" ■ (24) 



9 



It is interesting to note that Cartan's torsion tensor C^s = § ^§jpr vanishes when <7 M „ is Riemannian. Much of the 
mathematical literature has dealt with Finsler extensions of Riemannian geometry when the metric signature has 
been euclidean. However we need to consider pseudo-Riemannian structures. This can be formally done. However 
the Finsler norm leads to a metric which can lead to singularities in the metric for off-shell test particles. 

We focus on on-shell neutrinos (which are now known to have small masses). Again, this is motivated by our 
desire to discuss leptogenesis in such geometries. Moreover, for reasons that will become clear in section [V] it is 
neutrinos that play a preferential role in interacting non-trivially with the D-particle foam background, which induces 
stochastically fluctuating Finsler-like space times. As our main motivation is to compare the generic Finsler-like case 
with the D-foam model, as far as CPTV is concerned, we restrict our attention here on the effects of stochastically 
fluctuating Finsler geometries on dispersion relations of neutrinos and antineutrinos. We shall consider a particular 
type of Finsler metric on a manifold M which is known as the Randers metric [l4] 0. The norm Fr for the Randers 
metric is 

F R {x,y)=a(x,y) + p(x,y) (25) 

where 



a (x, y) = y/ (x) y^y v (26) 

and 

P(x,y) = b> l (x)y». (27) 

The conditions for a Finsler norm are satisfied. It was noted in [7l| that the geodesies of this metric coincided with 
the minimum time trajectories of a particle moving on a Riemannian manifold in the presence of a time independent 
drift given by a vector field. This is similar to Fermat's principle for propagation in refractive media. Similarities 
of D-particle foam to a refracting medium (oTI . I72I . [73| , will be mentioned briefly in the next section. If we were to 
assume that the result on minimum time trajectories was true for a pseudo-Riemannian situation and the drift was 
given by collisions due to D-particle scattering, then at a heuristic level a stochastic drift could be a reasonable generic 
phenomenological model of the back-reaction of low dimensional recoiling branes on matter. We shall write 

6 M (x) = (x) l„ (28) 

where i^is a constant vector. In our model <fi (x) will be a gaussian stochastic variable. On average the metric will be 
like a Riemannian metric if the mean of (<fi) vanishes. From (|2Tj) we deduce that 



, - -2 , -.7 1 i i r ^ v \ x ) r i- t Q \ X ) r vo \ x ) y e y° \ i / \ j c 

9V {x, y) = »V (x) + (f> (x) l^ly + r - -± -3 (x) l c y 

\a(x,y) a{x, y y J 

1 



-\ t r {t^q (x) y e (j) (x) l v + r VQ (x) y Q <p (x) l^) . 

a \x, y) 

We shall consider now a situation with (x) = r]^ where r/^ is the diagonal Minkowski matrix with entries 
(1, — 1, — 1, — 1). (The summation convention of repeated indices will be always understood unless explicitly stated 
otherwise.) Within the framework of a Robertson- Walker metric we shall ignore effects on the time-scale of the inverse 
expansion rate. Let us introduce some notation: 



a = vv^yV, 

and 



3 We mention, for completeness, that the other popular class of Finsler geometries, that ap pea rs in the General Relativistic version |54H57l | 
of the so-called Very Special Relativety Model [69], and cosmological extensions thereof [64, |64( are not characterised by CPTV in the 
dispersion relations, nor of the spin-curvature type discussed in section III CI In fact such VSR-related models have been propfosed in 
the past as candidates for the generation of Lepton- number conserving neutrino masses |70)|. and hence our Lepton- number violating 
considerations in this work do not apply. 
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We can then rewrite (x,y)as 



g»v = Vt*» + <f>%L + ( ^ - ) 0/3 + (C^ v + C u l M ) 05" 



However in order to construct energy-momentum dispersion relations we will need to calculate the dual metric h^ u 
(cf. (|24p) which involves y^ (u>) . We shall consider to be small and work to lowest order in the evaluation of G [x, uS). 
First we write y^ = y^ Q s + 0yf]\- From (l2"2")l we have 

w M = g^y (x, y( ) + 0y (1) ) (yfo + 0y ( l/ 1) ) . 
We can solve for 7/( )and y^.We find for j/( ) 

2/(0) = -"o, vlo) = ui, yf 0) = uj 2i yf Q} = w 3 . 

We find for y^ 

n _ l^w 2 - uj^d 



where w = y ti kS uj k ujs and <i = Z K u; K . Hence 

G(x,w) = n7-^/ K a; K (tn 2 + 2a;2). (29) 

We deduce that 

= rf v + cj)H^ (30) 

where 



for j = 1,2,3 



\2L 



2cj 2 ) + I k uj k (3m 2 (tu 



)-2wg(«7 2 -3w J a ))] > 



-ij (5ci7 4 wo + 2-n7 2 Wo) + (Zo ( ro4 + 2ro 2 u; 2 ) + l K u) K u)o (pvj 2 + 6a> 2 ))] 



and for i < j for z = 1,2,3 



H 11 = 



UJjh (tU + 2uj 2 



li (— m A — 6uj 2 uj 2 + vj 1 (2w 2 + w 2 )) (ro 2 — 6w 2 ) w l — ZjWj) 2lu 



117" 



We have assumed a homogeneous with being x independent. The mass shell condition Ww^Wv — m 2 leads to 
the following equation for luq: 

w2 ~ ~^ ^° ^ + + ll ^ ^ + Q ^ + 12 ^ 2 ^° + ^ + 13 ^ 3 ^ + ^ ^ = m2 ^ 

where up = lo\ + uj\ + ui 2 . 

In the model it is possible to choose / M . Not all choices will lead to asymmetric population distributions between 
particles and anti-particles. The space-like choice 1$ = gives a degenerate spectrum for particle and anti-particle and 
hence no CPTV in dispersion relations: 



w = 



20 
m 





1/2 r 


) +m 2 





1 (Ziwi + + Z3W3) 

m 



-1/2 



(32) 
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Therefore this case cannot be used for Leptogenesis in our framework. 
More generally the dispersion relation is 



oj = ±—l (2w 2 + m 2 ) +R(<f>,w,m) (33) 

TO 



where 



R (cj), uj, m) = (ti 2 + to 2 + ^ (2cj 2 + to 2 ) (huii + I2UJ2 + I3OJ3)) 



1/2 



The +sign is for the particle and the —sign is for the antiparticle. For the "time-like" case Iq = 1, l\ = I2 = h = the 
dispersion relation reduces to e 



CJ 



= yju> 2 + to 2 ± — (2u 2 + to 2 ) . (34) 



The sign of Iq can be reabsorbed in (j>. 

For the "null" case Iq = l\ = 1 and h — h — the dispersion relation reduces to 



w 



2 

m 



to V w' + m- 2 I I m 



y/& + m? l + £ 2- ^ 2 U ±^(2cJ 2 + to 2 ). (35) 



The parameters w play the role of momenta p in our case of neutrinos of mass to = to„ propagating in these space- 
times, and the Finsler metric may be seen as sort of back reaction on the space-time of such a propagation (to better 
appreciate this, the reader is invited to the discussion in the next section IVl where a particular model of D(efect)-foam 
is considered as a medium for neutrino propagation in the early universe, leading to Finsler-like metric distortions as 
a consequence of medum/particle interactions). 

Corresponding to such models involving D-foam, the parameter </> is modelled as a stochastic gaussian process with 
a mean a and standard deviation a. The interesting features in earlier analyses [21;], [22] were retained when a = 0. 
The fermion number distribution n from equilibrium statistical mechanics is given by 

f d 3 p 1 

n = 9d °t J (^expO^o-^ + l (36) 

where we have ignored degeneracy factors. In terms of spherical polars d 3 p = d^sinO dOdpp 2 where p = |pj , 6 lies in 
[0, 7r] and £ lies in [0,27r]. In the regime of validity of our analysis (i.e. relativistic neutrinos (pt ~ 0) that decouple 
at high temperatures T 3> to„ and for small stochastic fluctuations of the background geometry) we will consider 
ft = l/T (in units of the Boltzmann constant — 1) small and expand the relevant expressions in a power series. 
We denote n averaged over the distribution of <j) as <C n 3> . It is given by 



exp 



< n >= g d .o.f. / d(f> = n. (37) 

J-00 V^cr 

First of all, it immediately follows from (|3"4"|) , that for the "time-like" case, when a = 0, there is no particle/antiparticle 
asymmetry. This is to be expected, given that a 7^ corresponds in a sense to an averaged Lorentz violation in this 
stochastic geometry, and hence one of the basic assumptions for CPT Invariance of the effective theory of neutrinos 
in this "medium" is relaxed. 

For the "time-like" case, when a^O and /3 small, we obtain to leading order in T/m ^> 1: 

2 , / T \ f°° dxx 4 e x , 7tt 2 / T \ 

«An» agd . . f T 3 (-) =- agdo . f T 3 — - 38 

We require a < in order to have a particle-antiparticle asymmetry where the particle distribution dominates the 
antiparticle distribution. This yields the following Lepton (neutrino) asymmetry, assumed to freeze at the neutrino 
decoupling temperature ~ 10 15 GeV 

AL(T ~ T d ) = ~ —10 a — (39) 
s m v 

where, as usual, s denotes the entropy density, which for relativistic species is assumed to be s ~ gd.o.f.^- T 3 . 
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To obtain the phenomenologically correct value of AL(T ~ Tg) ~ 10 -10 , which is then communicated to the 
baryon sector via B+L violating sphaleron processes, or B-L conserving grand unified models (assumed appropriately 
embedded in such space-time geometries), one needs to take into account that, according to current data, the masses of 
the active neutrinos that are assumed to participate in (|39p must be smaller than m v < 0.2 eV. This implies then that 
one needs only an extremely small in magnitude violation of Lorentz symmetry on average in this stochastic Finsler 
space time, a ~ — 10~ 36 , in order to reproduce the observed Baryon Asymmetry in the universe. The assumption 
of fixing the sign of a is considered as fine tuning, and is a feature that is common in the models of gravitational 
leptogenesis/baryogenesis that exist in the current literature, as discussed briefly above [23l. \24L 126-32], 

Let us next calculate the asymmetry for the null case (|3l))) . For a = there is no asymmetry as we can see from 
considering (|37[) where we will exchange the order of the integrations over p and <f> . The (f> averaged expression for 
the particle distribution leads to 



where 



(6 — a) 2 \ T 3 P°° P 7 ' 1 
exp ( - vy ' — - / dxx 2 / dOsme- (1 + A(x,0)) (40) 



A(x,6) = { — [^x 2 {l + costf) + ^-z 4 (l + cos6) 2 } 

1 + e x e e- 



4</> 



2 



{l + e x ) 2 e 2 



x 4 {l + cos9) 2 } 



we have expressed p in terms of spherical polars, x = and the 9 — axis has been taken for convenience to 
coincide with the 1-axis; The corresponding expression for the antiparticle is 

d<t> exp ( - ^ ~ ^ \ / dxx 2 [ d9sin9— — (1 + B(x,9)) (41) 



<?Vk J-vo \ ° 2 ) (2tt) 2 Jo 



where 



B(x,0) = { [ — -x 2 (l - cos 9) + -rrx 4 (l - cos6») 2 ] 

1 + e x e e z 

4cj) 2 



e 

-x A {l-cos9) 2 } 



{l + e x f 

Although A(x, 9) ^ B(x, 9), for a = 0, on integrating over 9 (as part of the integration over p) we have 



< / d9sm9 A(x,9) > = < / d6 sin 9 B(x, 9) > 
Jo Jo 



Hence when a = there is no particle-antiparticle asymmetry. Hence we need Lorentz violation in the mean to obtain 
leptogenesis. This contrasts with our earlier work on correlations in neutral meson pairs created in meson factories 
where a signature of CPT violation [2l[ was present even for a = 0. 
When a ^ the counterpart of (|40]) is 



4 n x 



-P ( -^-3=- ) ( r 1 -;^ 2 ~ ~+—\ ) (42) 



(70? J_ 00 \ a 2 J \l + e x el 

and the counterpart of (|1T|) is 



ex P =5— h-r^^ +— t--;} • (43) 



a^/n J_ qo \ a 2 J \l + e x el 

Hence the analysis reduces to that for the time-like case. Consequently we have again (|38f and require a < to 
obtain an asymmetry with particle distribution exceeding that for the anti-particle. 

This need to tune the sign of the asymmetry in the approaches described so far leads us to pursue a different 
route to relations of the form (|12[) in the remainder of this paper . The route is not based on an approach involving 
effective local field theory. The framework we will follow fits naturally into a picture which has been advocated in 
the past to understand gravitational decoherence [2(| and dark matter abundance [74] . and may be considered to 
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be a more microscopic approach than the earlier proposals thatwe review. For D-foam the asymmetry is controlled 
by a parameter whose sign does not need adjustment. Furthermore, the constancy of the parameter, which has a 
stochastic origin, seems to be a more reasonable assumption than the same requirement for curvature. The D-particle 
foam model also has an interesting feature in that a non-Riemmanian metric is induced on the motion of matter 
on the brane world. Our mechanism of CPTV in D-foam did not directly rely on this aspect. However since this 
non-Riemannian structure is somewhat similar to the structure of Finsler metrics (see e.g. [51] and references therein) 
we have explored the implications for CPTV within the more general framework of stochastic Finsler metrics. We 
will find that it is indeed possible to have CPTV but in general it is necessary to choose the appropriate sign of s the 
mean of a stochastic parameter. 



V. STRINGY-DEFECT(D-)FOAM-INDUCED CPTV AND LEPTOGENESIS 



In this section we shall consider a population of DO-branes (or lower dimensional compactified D-branes which are 
effectively point-like from from the point of view of a brane world observer^) interacting with neutral fermions such as 
the neutrino and anti-neutrino. This interaction leads to different dispersion relations for neutrinos and anti-neutrinos 
which in turn leads to an excess of the population of neutrinos over anti-neutrinos. The freeze-out of neutrinos at 
the decoupling temperature of neutrinos leads to leptogenesis given by the standard cosmological considerations. The 
latter results, through standard Baryon (B) and Lepton (L) number violating sphaleron processes or B-L conserving 
interactions in grand unified models describing matter excitations on the brane, to the observed Baryon Asymmetry 
in the Unverse, with complete dominance of matter over antimatter, in a rather natural way, as we shall discuss below. 
Moreover, as we shall explain below, in this model of D(efect)-foam, the prevalence of matter over antimatter, i.e. 
the positive sign of the asymmetry An > 0, follows naturally, as a consequence of loss of energy of neutrinos during 
their interactions with the space-time defects, due to recoil of the latter. Thus, the sign of the induced asymmetry 
need not be fixed by hand, unlike the cases of gravitational leptogenesis discussed in previous sections. For instructive 
purposes, we first discuss the properties of the foam model, in the next subsection, before moving onto issues of CPTV 
and leptogenesis. 




FIG. 1. Schematic representation of a generic D-particle space-time foam model. The model of ref. fl6l. HH). which acts as a 
prototype of a D-foam, involves two stacks of D8-branes, each stack being attached to an orientifold plane. Owing to their 
special reflective properties, the latter provide a natural compactification of the bulk dimension. The bulk is punctured by 
DO-branes (D-particles), which are allowed in the type IA string theory. The presence of a D-brane is essential due to gauge flux 
conservation, since an isolated D-particle cannot exist. Open (F-)strings live on the brane world, representing Standard Model 
(SM) matter and they can interact in a topologically non-trivial way with the D-particle defects in the foam, only if they do 
not carry electric flux (electrically neutral excitations). Thus, from the SM matter excitations on the brane, mainly neutrinos 
interact non-trivially with the D-foam. Recoil of the D-particle during such interactions creates appropriate distortion in the 
space-time geometry, which depend on the momenta of the incident string states, and thus are of a generalized Finsler type. 
The propagation of neutrinos in such geometries is CPT Violating. 



4 It should be remarked that for the effective compactified D-"particles" the interactions with the charged matter excitations are suppressed 
relative to the neutral ones [73t . Hence, even in this case, it is the electrically neutral excitations which interact primarily with the 
D-foam. 



14 



A. The D-foam Model of the universe and Neutrinos 

D-foam models (l6l . ITsl l2l| are stringy models ol space-time foamy geometries, which involve brane universes, 
propagating in higher-dimensional bulk geometries. The bulk contains point-like D-brane defects ("D-particles" or DO 
branes) whose population density is constrained by the amount of CPTV that we observe. In many string theories 
(such as bosonic and type IIA string theories) they are stable zero-dimensional defects. However for our purposes 
we will consider them to be present in string theories of phenomenological interest [72| since, even when elementary 
D-particles cannot exist consistently, as is the case of type IIB string models, there can be effective D-particles formed 
by the compactification of higher dimensional D-branes [73[ (e.g. three-branes wrapped around three-cycles, with 
relatively small radii) . In general the construction of a model involves a number of parallel brane worlds with three 
large spatial dimensions, the required number being determined by target space SUSY.(Phenomenologically realistic 
models may require stacks of intersecting branes arranged in particular ways [25| .) These brane worlds move in a 
bulk space-time containing a gas of point-like bulk branes, called D-particles, which are stringy space-time solitonic 
defects [75| (c/. fig. [J). One of these branes is the observable universe. On this brane the D-particles will appear as 
space-time defects. Typically open strings interact with D-particles and satisfy Dirichlet boundary conditions when 
attached to them. Closed and open strings may be "cut" by D-particles, a process that involves capture of the incident 
open string and creation of stretched strings between the (recoiling) D-particle and the brane world (string "splitting") , 
and subsequent re-emission of the open string. It has also been speculated that nucleation of localised compactified 
defects |'2Gl | from such a D brane world (at the very high temperatures in the early universe) can be considered as a 
generation of compactified effective (metastable [76| but very long lived) D-particles. 

The preferential role of neutrinos in feeling the full effects of D-foam, and hence the CPTV, as we shall discuss 
below, is attributed to electric charge conservation: the representation of SM particles as open strings, with their ends 
attached to the brane worlds, prevents capture and splitting of open strings carrying electric fluxes by the D-particles. 
(We should recall that in string theory the electric charge is at the end point of an open string.) D-particles are 
electrically neutral and thus electric charge would not have been conserved if such processes had taken place. This 
is also consistent with the effective D-particles which may have formed as a result of nucleation [26[. Hence, the 
D-particle foam is transparent to charged excitations of the SM, leaving neutral particles, in particular neutrinos, 
susceptible to the foam effects. 

As discussed in detail in [22j the density of D-particles on the brane world is permitted to be relatively large, even 
at late eras of the universe, given the fact that bulk D-particles exert forces on the brane universe with mixed sign 
contributions to the brane vacuum energy, depending on the distance of the bulk D-particles from the brane (l7l . |771 | . 
Such forces are due to stretched strings between the defect and the brane. These energy contributions depend only on 
the transverse components of the relative velocities of the defect with respect to the brane worlds. In fact, depending 
on the distance of the bulk D-particle from the brane world, the sign of the contributions on the D-brane vacuum 
energy from the moving defect in the bulk, with velocity v perpendicular to the brane world, may be negative or 
positive. In particular, the interaction of a single D-particle, that lies far away from the D3 brane (D8-compactified) 
world, and moves adiabatically with a small velocity v± in a direction transverse to the brane, results in the following 
potential [3 

( longx 2 

V D0 -D8 = + 8to , » r » va'- (44) 

On the other hand, a D-particle close to the D3-brane (compactified D8), at a distance r' <C \/c7, moving adiabatically 
in the perpendicular direction with a velocity v^ hort will induce the following potential to it: 

hort W«^) 2 
VD0-D8 - ■ 

This difference in sign, then, implies that, one can arrange for the densities of far away and nearby bulk D-particles, 
which are not in general homogeneous, to be such that the total contribution to the brane world's vacuum energy is 
always subcritical, so that issues such as overclosure of the universe by a significant population of D-particle defects 
can be avoided. 

For our purposes in this work we may therefore consider that statistically significant populations of D-particles 
existed in the early eras of the brane universe. As the time elapses, the brane universe, which propagates in the 



5 For brevity, in what follows we ignore potential contributions induced by compactification of the D8 brane worlds to D3 branes, stating 
only the expressions for the induced potential on the uncompactified brane world as a result of a stretched string between the latter and 
the D-particle - the compactitication does not affect our arguments on the negative energy contributions to the brane vacuum energy. 
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higher-dimensional bulk (c/. fig. [I}, enters regions characterised by D-particle depletion, in such as way that the late 
eras cosmology of the universe is not affected. Nevertheless, as we shall discuss below, the early D-particle populations 
may still have important effects in generating neutrino-antineutrino populations differences (asymmetries) , which are 
then communicated to the baryon sector via the standard sphaleron processes [78] or B-L conserving GUT symmetries 
in unified particle physics models. 

To this end, we need to consider the effective dispersion relation of a (anti)neutrino field in a brane space-time 
punctured with statistically significant populations of D-particles. The latter is a dynamical population, consisting 
of defects crossing the brane all the time, thereby appearing to a brane observer as flashing "on" and "off" space-time 
"foamy" structures. The (anti)neutrino excitations are represented as matter open strings with their ends attached on 
the brane. The number density of (anti) neutrinos on the brane world is limited by the requirement that they do not 
overdose the universe, if neutrinos are assumed to have a chemical potential fi, then standard cosmological neutrino 
models predict that the number densities of a single flavour of relativistic neutrinos \emph{plus} antineutrinos in 
thermal equilibrium at temperature T v is estimated by [79] 

n3 3C 3 , 2ln2^ M 4 e. 



upon making the standard assumption that \i v <C T v for all neutrino flavours. The quantity £ v = is called 
the degeneracy parameter and is invariant under cosmic expansion, if we assume that the electron-neutrino chem- 
ical potential is the only one with significant presence in the early universe, then Big-Bang-Nucleosynthesis (BBN) 
constraints imply —0.04 < £„ c < 0.07. Thus, the order of magnitude of the neutrino plus antineutrino number 
density is agrees with naive standard estimate n v v ~ fx n 7i where n 1 is the photon density. Thus, today, where 
the temperature of the universe is of order To = 2.728~K (Cosmic Microwave Background), corresponding to an 
energy of kpTo ~ 2.35 x 10~ 13 GeV (with kp Boltzmann constant), the density of neutrinos is found to be of order 

— 112 cm -3 and scales roughly with the cubic power of the temperature: n yV ~ n\^l (j^J ■ So, for the decou- 
pling temperatures of neutrinos, kgT^ ~ 10 15 ~GeV, where we are interested in this work, in order to compute the 
frozen CPT Violating neutrino-antineutrino population differences, one obtains a number density of neutrino plus 
antineutrino populations of order 

n vV {T = T d ~ 10 15 GeV) - 10 85 cm' 3 . (47) 

On the other hand, as already mentioned, there are no similar restrictions on the population of the D-particle defects 
on the brane, in view of the negative contributions on the potential energy of the brane universe by bulk D-particle 
populations [22]. Thus, at the early universe, at the above neutrino-decoupling temperatures, we may even assume 
D-particle densities of one defect per string volume on the three brane world, without overdosing the universe. The 
assumption that the string length can take on values in the phenomenologically acceptable (post LHC era) range 
10 -27 — 10 -32 cm, corresponding to string mass scales from 10 TeV to 10 18 GeV, yields then a D-particle number 
density in the range 

n D (T = T d ~ 10 15 GeV) - 10 54 - 10 96 cm" 3 (48) 

respectively. Thus we observe that in order to be able to treat the D-particle populations as providing a more-or- 
less uniform "medium" over which neutrinos propagate, with non-trivial effective dispersion relations, we need to 
have at the decoupling temperature much higher densities of D-particles than those of neutrinos plus antineutrinos. 
Comparing (I47|) with (|48p . we observe that, if one assumes one D-particle per three-dimensional string volume on the 
brane, then this latter requirement excludes the low values of the string mass scale, implying an allowed range 

IQ-^Mp < M s < 10 _1 M P , (49) 

with Mp ~ 10 19 GeV the four-dimensional Planck mass. One of course could have much more dense D-particle gases 
in the early universe, which would allow for lower string scales. 



B. Kinematics of D-particle scattering and CPTV induced Leptogenesis 

We will now estimate the modification of the dispersion relations of neutrinos in such a "media" of D-particles in 
the early universe. The interaction of a string with a D-particle implies that at least one of the ends of the string 
is attached to the D-particle defect. Furthermore, the simultaneous creation of virtual strings stretched between the 
defect and the brane, describes the recoil of the D-particle. During the interaction time, the D-particle undergoes 
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motion characterized by non-trivial velocities, uu = 4j-Api — jj-TiPi along the brane longitudinal dimensions, where 
r, denotes the proportion of the incident neutrino momentum that corresponds to the momentum transfer Ap, during 
the scattering, and v± in directions transverse to the brane world [2(]| . 
As discussed in [111, H3, [z3l tne 

non-trivial capture and splitting of the open string during its interaction with the 
D-particle, and the recoil of the latter, result in a local effective metric distortion of the form: 

ds 2 = g^dx^dx" = (77^ + h^dx^dx" , h 0i = (u^^a a ) , (50) 

where u ; m is the recoil velocity of the D-particle on the D-brane world, with i = 1,2,3 a spatial space-time index, o~ a 
are the 2x2 Pauli flavour matrices with a = 1, 2, 3 (assuming two- flavour oscillations for simplicity). On average over 
a population of stochastically fluctuating D-particles including flavour changes, one may have the conditions (|55p . the 
second of which in the case of flavour oscillations can be generalised to 

< u l,i u l,j >= ^Sirfab ■ (51) 

(We still assume that <^ ujj , ~^>= . ) As a result of (fS"Tj) . on average, the flavour change during the interactions of 
neutrinos with the D-foam can be ignored. In such a case, any flavour structure in the metric (|5"0"|) is ignored 0. This 
result is the motivation for the consideration of a more general structure: Finsler geometry with stochastic parameters 
examined previously. 

However, the effects of D-foam go beyond those encoded in the induced Finsler like metric. The fine tuning that 
is required in stochastic Finsler metrics to get the correct sign for the particle-antiparticle asymmetry is a feature, 
although commonplace in other approaches, that is not entirely satisfactory. However because we have a microscopic 
model we can consider the kinematics of D-particle scattering. On considering string theory scattering amplitudes we 
find that the four momentum is conserved in the scattering of D-particles and strings. D-particles in the bulk exert 
forces on the vacuum energy of the brane world of mixed sign, depending on their relative distance. Thus, during 
the scattering process of a neutrino field with a D-particle, the vacuum energy of the brane fluctuates by an amount 
AV which depending on the process can be of either sign. From energy-momentum conservation, at each individual 
scattering event between a neutrino field and a recoiling D-particle, one could thus write: 

Pbcforc + Patter + — U|| = , ^before = Matter + ~ — u\ + AV (54) 

g s 2 g s 11 

where (p, -E)before (after) denote the incident (outgoing) neutrino momenta, energies repectively and we used the fact 
that the recoiling heavy D-particle of mass M s /g s (with M s the string scale and g s < 1 the string coupling, assumed 
weak, so that string perturbation theory applies) has a non-relativistic kinetic energy ^ ujj. We have also assumed 
that the fraction of the neutrino momentum transfer in the direction perpendicular to the brane world is negligible. 
The importance of the term AV not having a fixed sign in each individual scattering process is associated with the 
possibility of D-particle induced neutrino flavour oscillations [2l[ . 

Indeed, upon averaging ((...)) over a statistically significant number of events, due to multiple scatterings in a 
D-foam background, we may use the following stochastic hypothesis [21] 

<Cu i ||>=0, <C Uj || Uj ■ || <J 2 b~i 3 ■■ (55) 

implying that Lorentz invariance holds only as an average symmetry over large populations of D-particles in the 
foam. At a microscopic level, (|55p translates to momentum conservation on average in (|54|) . since -c u\\ ^>= 0. At 

an individual scattering process, if one represents the energy of the incident neutrino on-shell as ^Jp 2 + mf, where p 
is the amplitude of the conserved spatial momentum of the neutrino, and the outgoing one as yjji 2, + m\ , we observe 
that the energy-conservation equation (|54[) implies in general m% ^ mi. Which one is larger depends on the signature 



Ignoring the flavour structure, the metric 150t can be written as 

ds 2 = dt 2 + 2u l dx i dt - 8 ij dx i dx : > . (52) 
This metric was determined from world-sheet conformal field theory considerations [67l | and represents a dragging of the frame by the 
Gallilean (slowly moving) D-particle, which moves on a flat space-time background. However, the string excitations represent relativistic 
particles, and as such they move according to the rules of special relativity. Any four vectors attached to the strings, such as a four 
velocity, will evolve by a series of infinitesimal Lorentz boosts induced by the change of the D-particle velocity relative to the particle. 
In this sense, one may perform a time coordinate change in the metric l|52fl to write in in the form, up to terms u 3 for small recoil 
velocities |u| <g 1: 

ds 2 = dtff + 2u i dx i dt s - 5 ij (dx i - u i d%)(dx J ' - u j dt s ) + 0(u 3 ). (53) 
The metric 1531 is nothing but the so-called Gullstrand-Painleve metric [8Cj| , representing the geometry in the exterior of a Schwarzschild 
black hole, where the falling space into the black hole is represented as a Gallilean river on a flat space-time in which relativistic fishes 
swim. The river represents the frame of the recoiling D-particle, while the fishes are the relativistic matter strings. Here tg is the time of a 
free-floating observer who is at rest at infinity (compared to the centre of the black hole) . In the case of a black hole the relative velocities 
u l are coordinate dependent, of course, unlike our approximation in the D-foam case, although one may easily consider more general 
cases, where the recoil velocities of the D-particles in the foam are non uniform, in which case the analogy with the Gullstrand-Painleve 
river would become stronger. 
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of the term | u| + AV, which as mentioned is not of fixed sign, thereby allowing for neutrino oscillations to take 
place. The situation is somewhat analogous to the standard Mossbauer effect |8lJ, where the emitted or absorbed 
photon from a nucleus of an atom bound in a solid may sometimes be free of nuclear recoil, in contrast to the case of 
gases, thereby attributing the phenomena of nuclear resonances to such recoil-free fraction of nuclear events. In our 
case the r'i^cele of the "nuclei bound in a lattice" is played by the D-particle lattice. In addition to the D-particle recoil 
energy during scattering with stringy matter, which would lead to energy losses for the neutrinos, there are vacuum 
energy fluctuations, as a consequence of the motion of bulk particles in the foam, thus the neutrino experiences losses 
and gains from the vacuum, which results in the induced flavour oscillations. The analogue of resonances in this case 
would correspond to the loss-and-gain-free fraction of events, in which the neutrino does not oscillate. 

However, the effects of the D-foam go beyond the above-mentioned kinematical ones. On assuming isotropic 
momentum transfer, r, = r for all i = 1,2, 3. The dispersion relation of a neutrino of mass m propagating on such a 
deformed isotropic space-time, then, reads: 

PW = fa/** + V) = -™ 2 E 2 - 2E P- "II - f - ™ 2 = ■ (56) 
This on-shell condition implies that 

E = u» ■ p± \J (u ■ p) 2 + p 2 + m 2 . (57) 

We take the average <C • • • ^> over D-particle populations with the stochastic processes (|5ip , (|55p . Hence we arrive 
at the following expression for an average neutrino energy in the D-foam background: 



<-E> = <Cp-it>±« vV + m2 + (p-u) 2 > 

~ ±vV +m 2 (l + id 2 J , p>m, (58) 



for the active light neutrino species. The last relation in eq. (|58p expresses the corrections due to the space-time 
distortion of the stochastic foam to the free neutrino propagation. It is this expression for the neutrino energies 
that should be used in the averaged energy-momentum conservation equation (fM|) that characterises a scattering 
event between a neutrino and a D-particle. On further making the assumption for the brane vacuum energy that 
<C AV ^$>= 0, the total combined effect on the energy-momentum dispersion relations, from both capture/splitting 
and metric distortion, can then be represented as: 

^E 2 ^>=±^p 2 ~Tm 2 ~(l + )-A -IMl a 2 (59) 
V 2 / 2 9s 



Since antiparticles of spin 1/2 fermions can be viewed as "holes" with negative energies, we obtain from (|54|) and (|58|) 
the following dispersion relations between particles and antiparticles in this geometry (for Majorana neutrinos, the 
r'i^celes of particles /antiparticles are replaced by left /right handed fermions): 

1 ,\ 1M, 



2 \ _ __s Jl 

2 J 2 g s 



« £V » = vV + m 2 ( 1 + \a 2 ) + a 2 (60) 

V 2 / 2 gs 

where E > represents the positive energy of a physical antiparticle. In our analysis above we have made the 
symmetric assumption that the recoil- velocities fluctuation strengths are the same between particle and antiparticle 
sectors. (Scenarios for which this symmetry was not asssumed have also been considered in an early work (2l|.) There 
can thus be local CPTV in the sense that the effective dispersion relation between neutrinos and antineutrinos are 
different. This is a consequence of the local violation of Lorentz symmetry (LV), as a result of the non-trivial recoil 
velocities of the D-particle, leading to the LV space-time distortions (T5D|) . 

The discussion of CPTV in such foamy universes now follows the line of argument adopted by others: the difference 
in the dispersion relations between particles and antiparticles will imply differences in the relevant populations of 
neutrinos (n) and antineutrinos (n), (c/. the dispersion (|60l) ). This difference between neutrino and antincutrino 
phase-space distribution functions in D-foam backgrounds generates a matter-antimatter lepton asymmetry in the 
relevant densities 

« n - n »= g d . . f . J -0^ « [f(E) - f(E)} » , (61) 
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where Qd.o.f. denotes the number of degrees of freedom of relativistic neutrinos, and <C • • • 3> denotes an average over 
suitable populations of stochastically fluctuating D-particles (|55|) . 

Let us first make the plausible assumption that a 2 is constant i.e. independent of space. It is a parameter which 
can only be positive. Furthermore we will ssume that a 2 is independent of the (anti)neutrino energy. This is for 
estimation purposes only. We shall come back to a more detailed analysis later. We should note that the form of the 
dispersion relations (f5U|) is analogous to the case of CPTV axisymmetric geometries in the early universe, discussed 
previously (c/. (|12p ) with the ri^cele of the axial curvature scalar potential Bq being played here by the quantity 
i^cr 2 . In fact, it is easily seen that to leading order in a 2 the (1 + 5c 2 ) prefactors of the square roots in (|60|) play no 

r'i^cele, and hence the leading in a 2 contribution to the Leptonic asymmetry comes from the constant 5^ c 2 terms 
in the dispersion relations. The induced lepton asymmetry can therefore be calculated following similar steps as those 
leading to (|14p . upon the replacement of Bq by \^& 2 -, the difference being that here the value of Bq (i.e. of the 
D-foam recoil fluctuations a 2 ) is to be fixed phenomenologically. 

The result for the D- foam-induced lepton asymmetry can be estimated from (|6ip . using (|60p . Ignoring neutrino 
mass terms and (1 + ^-) square-root prefactors in (|60p . setting the (anti)neutrino chemical potential to zero (which 
is a sufficient approximation for relativistic light neutrino matter) and performing a change of variables \p\/T — > u we 
obtain from (f6Tj) the result: 

gd.o.f. „ 3 / ~2r 1 1 1 _ 9d.o.f. „ 3 



Au^^ff / duii 2 [ : : ] = £^£i T a Li 3 (-e-^)-Li 3 (- e ^ 

*^W^£)>0, (62) 



to leading order in a 2 , where in the last step we took into account the formal definition as a series of the Polylogarithm 

si 

--1 fc = 



function Li s (z) = X^fcLi W which is valid for \z\ < 1, while the cases \z\ > 1 are defined by analytic continuation. 



We thus observe that the CPTV term — |-^cr 2 in the dispersion relation (|6TJ)) for the neutrino, which corresponds to 
the energy 'loss' due to the D-particle recoil kinetic energies, comes with the right sign ('loss') so as to guarantee an 
excess of particles over antiparticles. Unlike the model of [24l I29I l3ll. l32| . then, where the sign of the Bo parameter 
had to be assumed, in our D-foam case there is no such freedom, and the positive An„ is derived from first principles. 
We consider this a nice feature of our model. 

As in standard scenarios of Leptogenesis, the Lepton asymmetry (|62p decreases with decreasing temperature up to 
a freeze-out point, which occurs at temperatures Td at which the Lepton-number violating processes decouple. This 
is taken to be the conventional one (in standard scenarios of Leptogenesis): Td ~ 10 15 GeV. 
The resulting lepton asymmetry then freezes out to a value (c/. Eq. (|14p ): 

AL(T < T d ) = — (63) 

s g s T d 

which survives today. The so calculated AL assumes the phenomenologically relevant order of magnitude of 10 -10 
provided 

^ a 2 ~ 10 5 GeV . (64) 

Since in these scenarios, the dimensionless stochastic variable, expressing fluctuations of a recoil velocity, is always 
less than one, a 2 < 1, one observes that the required Lepton asymmetry is obtained for D-particle masses larger than 

— > 100 TeV . (65) 
9s 

The so obtained AL can then be communicated to the baryon sector, t yield the observed (today) baryon asymmetry, 
via either B+L violating sphaleron processes, or B-L conserving interactions in Grand Unified theories, as in standard 
scenarios. 

The above estimate ignores an important fact, namely the dependence of the stochastic variable a 2 on the neutrino 
energy. Indeed, as discussed in detail in [20T [2l| . one may parameterise the momentum transfer by the fraction 
parameter of the incident momentum r, which is in turn assumed stochastic, that is 

Ui — -^-Apt — !• g s , no sum over i , <C r. t , <§; r,rj A 2 5ij . (66) 

1 V2 s Ms 
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In this case, the dispersion relations are modified by the replacement of 



2„2 



M? 



which is now momentum dependent: 



« E v » = ^Tm 2 ( 1 + ^A 2 p 



(67) 



« E? » = VP 2 + ( 1 + J^p A2 P 



2M? 
9 2 s 



9« 



2AL 

9s 
2M. 



■A 2 p 2 
■AV 



(68) 



Below, we shall evaluate the integral (|6 1 1) for the case (|55|) . assuming again A 2 <C 1 sufficiently small so that a 
truncation to order A 2 will be sufficient. A 2 will be assumed to be the same between particle and antiparticle sectors. 

In contrast to conventional point-like field theory models, where the upper limit of momentum integration can be 
extended to oo, in D-foam models, due to (|50p . this is extended up to the value for which the D-particle recoil velocity 
approaches the value of the speed of light in vacuo, c=l in our units, i.e. 



Pn 



\p\ max — 



9s %/A 2 

where r is the stochastic variable satisfying (|6"6")h The resulting integrals in (|6"Tj) then become: 



2lT Z 



dii. 



1 



1 



1 



- . - 2 g s A 2 T 



(69) 



(70) 



1 



The upper limit of the u integration cannot be taken as oo since An v should be evaluated at the decoupling temperature 
Td ~ 10 15 GeV, where it freezes out, as mentioned previously. An advantage of this CPTV as compared to the ones 
associated with axisymmetric geometries in the early universe [24], [29|, [3l], [HJ is that the CPTV parameter, A 2 , 
which characterises the fluctuations in the neutrino momentum transfer due to its interactions with the D-foam, is a 
parameter rather than a function. It is phenomenological (since it depends on the density of D-particles in the early 
universe) and cannot be significantly constrained by the cosmology of the early universe. We shall now analyse (|70D : 
first we set 



Tg s 



We note that 



An - 9d o f T 3 

_ 9d.o.f. T 3 

2tt z 



du 



I _|_ e u-u 2 A/2y 



du 



dv 



o l + e ^ _i£ 5- Jo 1 + e^+T- 



1 + e '«+« 2 A/2y 

dv 



(71) 
(72) 



We write T\ — J Q - 



dv 



+°xp(^-^f-) 



and I? 



Jo 



dv 



l+exp(^ + 



jja-. We will evaluate An„ in two limits:(i) a — > oo, A ~ 0(1) 



and (ii) A — > 0, — — > 0. For case (i) let x = 4r — \ and so 



dx 



(4 - v) 1 + cxp(ax) 



(73) 



where v is given by 



w(x) = 



1 - Vl - 2xA 2 



This integral can be evaluated (using Watson's lemma) to give 

, 1 A 2 3A 4 
Zi~A - + — + — + ... . 
a ar ar 



(74) 
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A similar analyis for I2 gives 



Hence Art„ in limit (i) is given by 



For case (ii) we have 



1 A 2 SA 4 

Zb~A(±-^ + ^- + ...). (75) 



IT 1 M, 



X x - f e~% v e^ and Z 2 - I e^e"^. (77) 
^0 Jo 

However the analysis again leads to the result (|7S)) . The lepton asymmetry resulting from ([75)1 freezes out at temper- 
ature Tj_ is 

AL(T<T d ) = ^ = ^^. (78) 
s M s 

We note that this result is compatible formally with if one takes into account (|67|) and associates the momentum 
maqnitude p of the realtivistic neutrino (i.e. its energy) with the temperature T . 

From ([75)1 , we observe that for a freeze-out temperature X^ ~ 10 15 GeV, the phenomenological value AL ~ 10~ ln 
is attained for 

^~ 10 25 A 2 GeV. (79) 
9s 

For A 2 ~ 10~ 6 a Planck size D-particle mass M s /g s ~ 10 19 GeV is required so that the D-foam provides the 
physically observed Lepton and, thus, Baryon Asymmetry. For the unnaturally small A 2 < 10~ 21 one arrives at 
M s /g s ~ 10 TeV. Unfortunately, for A 2 ~ O(l) transplanckian D-particle masses are required. We should stress 
that the above conclusions were based on the assumption that the freeze-out temperature was the temperature at 
decoupling of neutrinos in standard big-bang cosmology. 

Our approach to leptogenesis is distinguished from others in that a local effective field theoretical description is 
not adopted. Because of D-particle recoil when scattering off matter strings the background of D-particles can be 
modelled as a stochastic medium [20T-[22j. The underlying string theoretic description provides the rigorous description 
of the scattering of D-particles. The D-particles backreact (as seen from infra-red divergences in perturbation theory) 
and change the metric which influences the space in which matter is moving. Furthermore as discussed at length in 
[22l |. and mentioned briefly above, the D-particle foam model does not lead to overdosing the universe. Hence despite 
having statistically significant populations of D-particles in the early universe, which provide the CPTV background 
on which neutrinos propagate, the assumption of a subcritical energy density for the universe can still hold. 



VI. CONCLUSIONS AND OUTLOOK 



In this work we have considered some models leading to CPTV gravitational leptogenesis, based on string-inspired 
constructions. The primary model that we discuss here involves a brane universe propagating in a higher-dimensional 
bulk space-time, which contains populations of ("point-like" ) DO brane defects (D-particles). The topologically non- 
trivial interaction of the D-particles with the neutrinos, in the sense of capture and subsequent re-emergence of the 
latter by the defects, is allowed by electric charge conservation; in contrast the medium of defects is "transparent" for 
charged excitations of the standard model (such as charged leptons and quarks). The propagation of neutrinos in this 
medium of space-time defects results in CPTV dispersion relations, which are different for neutrinos and antineutrinos. 
It is essential to notice that this difference is obtained by considering the kinematics of the defect /neutrino scattering 
on the assumption that the antiparticle has negative energies following the"hole theory" of Dirac. These considerations 
apply only to fermions, as opposed to bosons, because the "hole theory" is based on the fermion exclusion principle. 
This model has two aspects. One is the induced metric in space-time, which depends on momenta as well as space- 
time coordinates, and the other is the particle-like kinematics involved in the DO-brane scattering of matter. The first 
aspect is similar to stochastic Finsler metrics which we have considered in a general way without necessarily tying 
it to the brane model. The second, kinematical aspect is important in leading naturally to matter dominance over 
antimatter, in contrast to earlier proposals of gravitational leptogenesis, where this sign has to be adjusted by hand. 
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The cosmology and particle-physics phenomenology of such scenarios , although currently in their infancy, are worth 
pursuing in our opinion. In particular, for the defect model, at early epochs of the universe, despite the fact that 
significant populations of massive D-particles are assumed present, we have argued here and in our previous works 
[22I |74| that one can avoid overclosure of the brane universe: the defects are allowed to propagate in the bulk and 
thus they exert forces on the brane world which are such that there are mixed sign contributions to the brane vacuum 
energy, depending on the distance of the bulk defect populations from the brane. The population of D-particles 
in the bulk is also constrained by the amount of Lorentz violation at late eras of the universe, which leave their 
imprint on the Cosmic Microwave Background [82, 83] or in vacuum- refractive- index tests of arrival times of cosmic 
photons [H,[li-ii|. 
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